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Abstract 

Ground-state phase diagram of the mixed spin- 1/2 and spin-1 Ising chain with axial and rhombic 
zero-field splitting parameters is exactly calculated within the framework of the transfer-matrix 
method. It is shown that the rhombic zero-field splitting parameter prefers the magnetically ordered 
phase instead of the disordered phase. 
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I. INTRODUCTION 



Exactly solved one-dimensional quantum spin models traditionally belong to the most 
fascinating research they provide valuable insight into otherwise hardly under- 

standable aspects of cooperative and quantum phenomena jl| . In this work, we will exactly 
treat the mixed spin- 1/2 and spin-1 Ising chain with axial and rhombic zero- field splitting 
parameters. 

II. MODEL AND ITS EXACT SOLUTION 

Consider the Ising model for a chain consisting of the alternating spin-1/2 and spin-1 
atoms, which accounts also for axial and rhombic zero-field splitting parameters. The total 
Hamiltonian of this spin system can be written as a sum of two terms i-L = T-Lex + 'H^zfs- The 
former term accounts for the Ising-type exchange interaction between the nearest-neighbor 
spins 

N 

iiex = -JY.si{^l + K+^l (1) 

fc=l 

and the latter term accounts for the axial {D) and rhombic [E) zero-field splitting (ZFS) 
parameters acting on the spin-1 atoms only 

N N 

^i/l = -D Y.iSlf - EY.[{Slf - {Sin (2) 

fc=l k=l 

Above, al and S^{oi = x,y,z) denote standard spatial components of the spin-1/2 and 
spin-1 operators, respectively, denotes a total number of spin-1/2 (spin-1) atoms and the 
periodic boundary condition a^^i = ai is imposed for further convenience. It is worthwhile 
to remark that there is one-to-one correspondence between the Hamiltonian 'H^zfs given by 
Eq. ([2]) and the Hamiltonian with three different parameters D^, and 

N N N 

= E(^D' - T.iSir - D^Y^iSlf. (3) 

k=l k=l k=l 

The equivalence between ii^^fs ^"^^ '^ffs can easily be verified by establishing a rigorous 
mapping correspondence between the relevant interaction terms appearing in the Hamilto- 
nians ([2]) and ([3]). The total angular momentum of the spin-1 atoms is integral of motion 
= (-S*^)^ + {Siy + {S^y = 2 and hence, one of three parameters D^, and must 
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depend on the other two. Consequently, the Hamiltonians i-i^zfs ^^"^ ^i/s differ one from 
the other just by some constant factor ii^^fs — ^^zfs + ^i whereas the relevant interaction 
terms C, D and E are connected to the ones , and through the mapping relations 

C = D- + Dy,D = D^- , E = ■ (4) 

The model under investigation thus turns out to be equivalent to the one recently studied 
by Wu et al. [21] using the approach based on Jordan- Wigner transformation. 

Here, the investigated model system will be exactly treated within the framework of 
transfer-matrix method First, it is useful to rewrite the total Hamiltonian as a sum of 
site Hamiltonians ik. = J2k^k, whereas each site Hamiltonian involves all the interaction 
terms associated with the spin-1 atom from the kth lattice site 

= -JSl{al + al^,) - D{Slf - E[{Slf-{Sin (5) 

Due to a validity of commutation relation between different site Hamiltonians, the partition 
function can be partially factorized into the product 

N 

^ = E nTr5.exp(-/3Hfc), (6) 

{ak}k=l 

where /3 = l/{k-BT), ks is Boltzmann's constant, T is the absolute temperature, Tr^^, means 
a trace over spin degrees of freedom of the kth spin-1 atom and Y^{ak} denotes a summation 
over all possible configurations of the spin-1/2 atoms. After tracing out spin degrees of 
freedom of the spin-1 atom, the relevant expression on r.h.s of Eq. will depend just on 
its two nearest-neighbor spins ak and ak+i- Moreover, this expression can be subsequently 
used in order to define the transfer matrix 

Jlak, (Tfc+i) = TiSk exp{-f3'Hk) 

=1 + 2 exp(/3D) cosh [l3p\al + a^+J^ + _ (7) 

The rest of our exact calculations can be accomplished using the standard procedure devel- 
oped within the transfer- matrix approach [3']. This rigorous technique allows one to express 
the partition function in terms of respective eigenvalues of the transfer matrix 

N 

^ = E n = TrT^ = + A''- (8) 



In the thermodynamic hmit N ^ oo, the free energy per unit cell can be expressed solely 
in terms of the largest eigenvalue of the transfer matrix 



where Tn = T(±l/2,±l/2) and T12 = T(±l/2, =f1/2) were used to denote two different 
matrix elements of the transfer matrix defined through Eq. (^^. 

III. RESULTS AND DISCUSSION 

Now, let us take a closer look at the ground-state behavior of the investigated model 
system. For simplicity, our subsequent analysis will be restricted only to the particular case 
with the ferromagnetic interaction J > 0, since the relevant change in sign of the parameter 
J causes just a rather trivial reversal of all the spin-1/2 atoms. 

In the zero temperature limit, the first-order phase transition line given by the condition 
D = —\\P + E^ separates the ferromagnetically ordered phase (OP) from the disordered 
phase (DP). The relevant spin order appearing in the OP and DP can be unambiguously 
defined through the eigenvectors 



where the product runs over all lattice sites, the former (latter) ket vectors specify the 
state of the spin-1/2 (spin-1) atoms and the mixing angle is given by = arctan(ii^/ J). 
In the DP, all the spin-1 atoms tend toward their 'non-magnetic' spin state |0) on behalf 
of a sufficiently strong (negative) axial ZFS parameter and hence, each spin-1/2 atom may 
completely independently choose any of two available spin states |±l/2). However, the more 
striking spin order emerges in the OP, where the magnetic behavior of the spin-1 atoms is 
governed by a quantum entanglement of two magnetic spin states |+1) and |— 1) and all 
the spin-1/2 atoms reside their "up" spin state |l/2). In this respect, the rhombic ZFS 
parameter gradually destroys a perfect ferromagnetic order between the spin-1/2 and spin-1 
atoms, which appears in an absence of the rhombic term. 

For better illustration, Fig.[T]^a) depicts the ground-state phase diagram in the E/ J—D/ J 
plane. The most surprising finding stemming from Fig. [I](a) is that the phase boundary 




kBT\n{Tn+Tu), 



(9) 
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FIG. 1: Ground-state phase diagram in two different planes: (a) E/ J — D/ J plane; (6) 2J/Dx — 
Dz/ J plane for Dy = 0. 

between OP and DP shifts toward more negative values of the axial ZFS parameter when 
increasing a strength of the rhombic ZFS parameter. Accordingly, it turns out that the 
quantum entanglement between the spin states |+1) and |— 1), which is caused solely by 
the rhombic ZFS parameter, energetically stabilizes the OP before the DP. For comparison. 
Fig. [U^b) illustrates the ground-state phase diagram in the U jD^ — D^j J plane when using 
Eq. ([3]) in order to define the ZFS Hamiltonian. Note that this phase diagram is in accord 

I i 

with the recent results of Wu et al. [2], but this phase diagram is apparently less convenient 
for interpreting the phase boundary between OP and DP as the parameter changes 
according Eq. both axial as well as rhombic ZFS parameters. 

In conclusion, it is worthy to notice that the rigorous procedure developed on the grounds 
of the transfer-matrix method can readily be adapted to treat the investigated model system 
even in a presence of non-zero external magnetic field, which will be examined in detail in 
our forthcoming work. 
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